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ABSTRACT 

Variations of the microturbulent velocity with phase and 
height in the atmosphere have been reported in classical Cepheids. 
It is shown that these effects can be understood in terms of 
variations of the velocity gradient in the atmospheres of these 
stars. 
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I. INTRODUCTION 

Struve (1932) introduced the concept of microturbulence to 
explain the anomalously large Doppler broadening velocities he 
found for supergiants. This was explained by Struve and Elvey 
(1934) as being due to either a "turbulence of small eddies" or 
to "several shells which expand with different velocities." It 
is now known that microturbulence varies with height in the 
atmospheres of supergiants (Wright 1946, Huang and Struve 1960) 
and with phase in Cepheids (van Paradijs 1972). Differential 
motions have been reported in the atmospheres of supergiants of 
many spectral types (Abt 1957, Ayden 1972, Rosendhal and Wegner 1970) 
and in particular in Cepheids (van Hoof and Deurinck 1951, Dawe 1969). 
Van Paradijs (1972) has noted that the microturbulence in Cepheids 
is a maximum near the phase of most rapid contraction. Dawe (1969) 
and van Hoof and Deurinck (1951) have shown that the velocity gradient 
is appreciable at this phase. The problem of the effect of a velocity 
gradient on the curve of growth was investigated by Kubiowski and Ciurla 
(1965) and Ciurla (1966) for a B2V star. Their work is extended in 
this paper to a cooler supergiant model. It will be shown that the 
microturbulence can be explained by a velocity gradient in the 
atmosphere that is consistent with observed velocities of the 
spectral lines. 
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II. METHOD 

The method for computing line profiles in a moving atmosphere 
follows a suggestion made by Chandrasekhar (1945). The line 
absorption is described by a Voigt profile, H(a,u), where a_ is the 
damping parameter arid u = (A-A 0 )/AA D . When there are velocities 
in the atmosphere, u_ must be modified to account for the motions 
of the gas. This can be done by letting u' = u + yvA/cAA D f where yv 
is the local velocity of the gas projected onto the line of sight, 
and using H(a,u') to compute the line opacity. The specific intensity 
can be computed from 

I v (0,y) = /\ (t v ) e"V y dt v /y 
and then the flux from 

F (0 ) = 2 I'll 0,y)ydy. 

V O V 

Other methods commonly used, such as the Feautrier method or the 
quadrature integration of Milne's second equation (Kourganoff 1952), 
cannot be used since the line opacity is a function of y . 



III. RESULTS 


To make the test case as realistic as possibles a model 
atmosphere from Parsons (1969) and an Fel line were used. The 
model has T ^ = 6300°K and log g = 1.8 while the 4494. 57$ Fel 

line (excitation potential 2.2 ev) was chosen to give a reasonable 

variation of the number of absorbers with depth in the atmosphere. 

To study a wider range of effects than Ciurla (1966), curves 
of growth were computed for log a = -1, -2, -3 and microturbulent 
velocities £ = 0 and 5 km s”"* by varying the number of absorbers 
in the line of sight. These curves are given in Figures 1-3 and 

table 1. (The curves for £ = 0 and £ . = 5 km s~^ do not come 

together at large n/n_ because the ordinate is -log W/A instead of 
“log W/AAp.) Underhill (1947) has shown that a velocity of expansion 
(or contraction) constant in t cannot change the equivalent width, W, 
of a line. Such a velocity field will produce asymmetric line 
profiles, however, due to the integration over the surface. As 
a check several profiles were computed with v (t) = 20 and 40 km s - ^.. 
In no case was the change in W greater than 1%. This change is due 
to errors in the angle and frequency integration and can be used as 
a crude estimate of the errors in all these calculations. 

Since little is known about the velocity field in the atmosphere 
of a Cepheid, an arbitrary choice of_v (t) was made, v (t) = -a log t. 
This is convenient because it allows a correlation of mean optical 
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depth of formation of a line and its observed radial velocity 


and is nearly linear with geometrical height in the line forming 
region. An arbitrary constant may be added to y^(x) without 

changing W, but it will change the shape of the profile. 

The results for 5 = 0 and a = 5 and 10 km s~^ are shown in 

Figures 1-3 and Table 2. With a = 10 km s”^ the curve of growth 

( Cctl 0 ) ls near ty identical to the normal curve of growth with 5=5 km s - ^ 
(C^g) until the damping portion is reached. In all cases C^g has 
a wider plateau than C^. An observer would interpret this as 
being due to a lower value of the damping parameter a^ The decrease 
in a^ at phases when 5 is large has been observed by Rodgers and 
Bell (1968a, 1968b). This change in a^ is easily understood. The 


vertical shift between the damping parts of and C^g is propro- 
tional to the ratio of the Doppler widths. Since C^g and C^g 


have the same Doppler width as C^g, in the strong line asymptotic 
limit the three curves must join. The only way this can happen 


is for C^iq to be below C^g as the lines get strong. 

Ciurla (1966) has done a similar calculation for a B2V model 
and a line with log a^-1. His results agree well with those 
presented here but, he needed a larger velocity gradient to mimic 
5 = 5 km s~^. The velocity difference between the strongest and 
the weakest lines is about 17 km s' 1 as opposed to the 9 km s“^ 


obtained here. 
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Dawe (1969) has plotted observed velocity versus mean optical 
depth of formation for weak lines in l Carinae. Reading from Dawe's 
figure 3, the weakest lines, formed near t = 0.3, show a velocity of 
about 18 km s"^ while those formed near x =0.1 show about 22 km s . 

I use a = 15 km s~^ to correspond roughly to £ = 7.5 km s”^ observed 
by Rodgers and Bell (1968a) near this phase. Using v = -alogx, 
v(0.1) - v(0.3) = 7.2 km s - ^ . Correcting for the integration over 
the surface by the factor v i $ / v rac j = 24/17, the predicted velocity 
difference is 5.0 km s”"*. The radial velocities "observed" from 
the minima of the computed profiles are given in Table 3. An exact 
comparison is not meaningful since the observed velocities and shapes 
of the lines are more sensitive to the velocity distribution than 
is the shape of the curve growth. An attempt to match the 
observed velocities and profiles more exactly must wait until more 
is known of the true velocity field. 
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IV. CONCLUSIONS 

It has been shown that at least part of the microturbulence 
in Cepheids can be explained by a velocity gradient in the 
atmosphere that is consistent with the observations. Variations 
in the velocity gradient with phase and height in the atmosphere 
can produce the observed variations in microturbulence. 

The method can also be applied to normal supergiants to see 
if their microturbulent velocities can be explained in the same 
way. Rosendhal (1970) has shown that microturbulence decreases 
from AO to later spectral types in supergiants and Rosendhal and 
Wegner (1970) have reported velocities in the atmospheres of A 
supergiants which they associate with mass loss. If this inter- 
pretation is correct, the Oort constant A and the possible expansion 
term K determined from radial velocities of supergiants will have a 
systematic error. Due to the expansion of the atmosphere, all the 
measured radial velocities will be too small by several km s”^ . 
Discrepancies between optical and radio determinations of the spiral 
structure of the galaxy have been discussed and are referenced by 
Mihalas and Routly (1968). If the variation in microturbulence 
with spectral type can be interpreted as a variation in the velocity 
of expansion of the atmosphere, some of the scatter in the determination 
of the galactic structure constants can be explained. 
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TABLE 1 


-log W/A vs log of Relative Number of Absorbers 
No Velocity Gradient 


log a , 
£(km s' 1 ) 
log n/n 0 

-1 

0 

-1 

5 

-2 

0 

-2 

5 

-3 

0 

-3 

5 

-1.000 

5.203 

5.023 

5.235 

5.031 

5.241 

5.031 

-0.699 

5.060 

4.795 

5.110 

4.811 

5.117 

4.812 

-0.399 

4.949 

4.610 

5.022 

4.637 

5.032 

4.641 

-0.097 

4.850 

4.470 

4.953 

4.513 

4.968 

4.518 

0.204 

4.747 

4.356 

4.894 

4.423 

4.917 

4.432 

0.505 

4.634 

4.257 

4.833 

4.354 

4.874 

4.368 

0.806 

4.511 

4.154 

4.776 

4.295 

4.836 

4.317 

1.107 

4.373 

4.034 

4.693 

4.238 

4.802 

4.274 

1.408 

4.234 

3.922 

4.616 

4.179 

4.768 

4.237 

1.709 

4.093 

3.797 

4.510 

4.115 

4.727 

4.203 

2.010 

3.953 

3.670 

4.398 

4.035 

4.683 

4.171 

2.311 

- 

- 

- 

- 

4.625 

4.138 

2.612 

- 

- 

- 

- 

4.536 

4.102 

2.913 

- 

- 

- 

- 

4.415 

3.944 
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TABLE 2 

-log W/A vs log of Relative Number of Absorbers 
v(t) =-a log t 


log a 
a(km s" 1 ) 
log n/no 

-1 

5 

-1 

10 

-2 

5 

-2 

10 

-3 

5 

-3 

10 


- 1.000 

5.138 

4.994 

5.099 

5.039 

5.091 

5.040 


0.204 

4.528 

4.392 

4.608 

4.429 

4.621 

4.432 


1.408 

4.174 

4.083 

4.417 

4.222 

4.454 

4.240 


2.612 

3.666 

3.635 

4.119 

4.027 

4.334 

4.147 




Radial 

Veloci ty 
Profiles 

TABLE 3 

in km s~^ Observed from Computed 
for Cases in Table 2 

log a 

-1 

-1 

-2 -2 

-3 

-3 

a(km s' 1 ) 

5 

10 

5 10 

5 

10 

log n/n 0 






-1.000 

2.9 

5.4 

3.3 5.4 

2.1 

5.4 

0.204 

5.4 

6.6 

5.4 9.3 

4.2 

9.3 

1.408 

6.7 

13.3 

6.6 14.7 

6.6 

14.7 

2.612 

10.8 

14.7 

9.3 14.7 

6.6 

14.7 
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CAPTIONS TO FIGURES 


Fig. 1. 


Fig. 2. 
Fig. 3. 


Curves of growth for log a = -1. — “ . — , 4=0 (C^ Q ); 


» ^ — 5 km s ( C^. ^ ) 9 o~ 

■i 


-o, a = 5 km s 


-1 


(C a5 ); • #, a = 10 km s (C^q). ^ t * ie m '' cro " 

turbulent velocity and a is the velocity gradient parameter. 
Curves of growth for log a = -2. Notation is the same as 
in Figure 1. 

Curves of growth for log a = -3. Notation is the same as 
in Figure 1. 
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